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Abstract 

u i 

We address the optimization of the sum rate performance in multicell interference-limited single- 
hop networks where access points are allowed to cooperate in terms of joint resource allocation. The 
resource allocation policies considered here combine power control and user scheduling. Although very 
promising from a conceptual point of view, the optimization of the sum of per-link rates hinges, in 

(N 

transmitter channel information feedback across all network cells. In this paper, we show that, in fact, 
distributed algorithms are actually obtainable in the asymptotic regime where the numbers of users per 



principle, on tough issues such as computational complexity and the requirement for heavy receiver-to- 



cell is allowed to grow large. Additionally, using extreme value theory, we provide scaling laws for 
upper and lower bounds for the network capacity (sum of single user rates over all cells), corresponding 
to zero-interference and worst-case interference scenarios. We show that the scaling is either dominated 
by path loss statistics or by small-scale fading, depending on the regime and user location scenario. 

We show that upper and lower rate bounds behave in fact identically, asymptotically. This remarkable 
result suggests not only that distributed resource allocation is practically possible but also that the impact 
of multicell interference on the capacity (in terms of scaling) actually vanishes asymptotically. 



Part of this work was presented at the 3rd IEEE workshop on Resource Allocation in Wireless Networks (RAWNET07), 
Limassol, Cyprus, April 16th 2007. 
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I. Introduction 

The performance of wireless cellular networks is impaired by the problem of interference. 
Traditional ways to tackle this problem include careful planning of the spectral resource and the 
use of interference mitigation or advanced coding/detection techniques combined with fast link 
adaptation protocols at the physical layer [1], [2]. In a typical approach to resource planning, the 
system designer aims at the fragmentation of the network geographical area into smaller zones 
(reuse patterns), which are rendered more or less isolated from each other from a radio point 
of view by the assignment of slices of the spectral pie which are mutually orthogonal. In this 
case, interference is reduced, but at the expense of some of the overall spectral efficiency. More 
subtle resource allocation protocols include power control [3] and dynamic channel assignment 
methods which also help alleviate the problem of interference. However the majority of such 
techniques are initially rooted in voice-centric network design and, as such, are geared toward 
achieving a given signal to noise plus interference ratio (SINR), common to all users, rather 
than maximizing the spectral efficiency defined as the number of Bits/Sec/Hz per area [4], [5]. 

As an alternative to such techniques, a recently proposed approach to enhancing the network 
overall throughput consists in turning interfering links into active links capable of carrying infor- 
mation to the user. In such as scheme the transmitters cooperate to offer a virtual MIMO setting 
[6], [7], [8]. Although cooperation involving joint processing and encoding at all interfering 
transmitters constitute an efficient weapon to enhancing the system capacity, it also presents the 
designer with practical challenges such as limiting the overhead of cell-to-cell communication 
and the need for symbol- accurate synchronization between the transmitters. 

Another important limitation common to most of the above mentioned work is the fact 
that interference mitigation and joint transmitter processing are addressed, as physical layer 
enhancements, while ignoring the possible impact of the user scheduling algorithm implemented 
in the multiple access protocol. We already know that single-cell greedy user scheduling has a 
profound impact on the fading statistics seen by the above layers of the system. In the case of 
multicell user scheduling, we may naturally wonder if interference may be impacted as well. In 
other words, a broad question is what is the impact of interference when a scheduling protocol 
is available that allows us to select within a pool of users? 

Thus in this paper, we revisit the problem of interference in single-hop cellular wireless 
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networks and examine it in the light of multicell resource allocation algorithms. The considered 
algorithms involve some level of cooperation between the transmitters in the form of joint power 
control and user scheduling, with the aim of maximizing of the network throughput. Since joint 
processing/encoding is not considered here, in a first approach, the network throughput is defined 
in an information theoretic sense as the sum over all cells of the per-cell single-user decoding 
rates, achieved by simultaneously receiving users. 

A simple intuitive idea behind multicell resource allocation is to exploit the large amount of 
spatial and multiuser diversity offered by the extra multicell dimension in order to optimize the 
network performance at all times. Clearly the potential gains comes with great challenges. One 
is the complexity associated with the joint optimization of a large number of parameters (slot 
assignment, power levels, ..). Another one is the need for the joint processing of multicell channel 
state information which necessitates huge cell-to-cell signaling overhead. This makes global 
network coordination hard to realize, especially in fast mobile settings. Despite the challenges, 
some recently published work suggests possible techniques for low complexity and distributed 
resource allocation. Examples of such approaches include game theoretic algorithms with pricing 
[9], [10], resource allocation based on quantized power levels, and iterative/greedy capacity 
maximization techniques. An overview of such techniques is available in [11]. Nonetheless it 
remains that such approaches are suboptimal and do not help characterize the system behavior 
when the number of user is large. 

In this paper we address the problem of sum rate-maximizing resource allocation, where 
the protocol takes the form of joint multicell power control and scheduling. We investigate 
upper and lower bounds on the maximum network capacity provided by resource allocation in 
interference-free and full-powered interference scenarios, respectively. Interestingly the solution 
to the multicell scheduling and power control is fully distributed in both scenarios. We study 
these bounds in the asymptotic regime where the number of users per cell is allowed to grow 
large while the number of cells remain fixed. We introduce scaling laws of capacity for this 
asymptotic regime, based on extreme value theory. We show that for certain idealized network 
settings, our models coincide mathematically with recently published models in the different 
context of opportunistic single cell beamforming, in particular those of e.g. [12], [13]. Therefore 
certain results from this domain can be reused. For other general network settings, different 
results of extreme value theory must be explored. 
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After identifying the scaling laws for the upper and lower bounds on capacity for two different 
forms of network/channel models, we draw the surprising conclusion that these are in fact 
identical, in both cases. In the particular case of users located at random in a finite disk around 
each access point, the analysis reveals a much higher growth rate than in the case of users having 
an equal average SNR. 

One of the practical impacts of our results is the suggestion that very simple distributed 
resource allocation algorithms can be used with little loss of performance on the network capacity 
compared with a centralized optimal scheme. Another interesting lesson is a better understanding 
of the the price paid due to multicell interference in terms of rates, when the number of users 
becomes large and when a maximum-rate scheduler is exploited. Our theoretical claims are 
backed with Monte-Carlo simulations on multicell fading channels. 

II. Network and signal models 

We consider a wireless network featuring a number of transmitters and receivers. Among these, 
there are N transmit-receive active pairs, which are simultaneously selected for transmission by 
the scheduling protocol at any considered instant of time, others remaining silent. All active links 
interfere with each other. This setup, an instance of the interference channel [14] can be observed 
in e.g. a cellular network with reuse factor one, such as the upcoming 802.16 (WiMax) and 3GPP 
(LTE) wireless standards. We assume each of the iV cells is equipped with an access point (AP) 
and that APs communicate with the users in a single-hop fashion. We also assume the APs are 
time-synchronized. In this paper we focus on the performance of downlink communication from 
the AP to the users. However we believe our analysis carries over to the uplink without great 
difficulty. 

Let U n be the number of users randomly distributed over cell n, for n = 1, . . . , N. We will 
assume these users are uniformly randomly distributed over either a circle or a disk around their 
access point. 

Within each cell, we consider an orthogonal multiple access scheme so that a single user is 
supported on any given spectral resource slot. A resource slot can be a time or frequency slot in 
TDMA/FDMA, or a code in orthogonal CDMA. For instance in OFDMA-based WiMax or LTE 
standards, a resource slot is represented by a unique time/frequency slice. For ease of exposition, 
single antenna devices are considered. Generalizations to MIMO-aided spatial division multiple 
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access are currently under investigation. On any given spectral resource slot, shared by all N 
cells, we denote by u n G {1, . . . , U n } the index of the user that is granted access to the slot (i.e. 
scheduled) in cell n. An example of such a situation is depicted for a simple two cell network 
in FigfTJ 

We denote the complex downlink channel gain between the z-th AP and user u n of cell n 
by a Unj i. The local channel state information (CSI) is assumed perfect at the receiver side. This 
information is also fedback to the control unit responsible for resource allocation, either in a 
centralized or distributed manner (this point crucial when it comes to applicability, as discussed 
later). The received signal Y Un at user u n is given by 



where X Un is the message-carrying signal from the serving AP, subject to a peak power constraint 
Pmax- J2i^ n a u n ,iX Ui is the sum of interfering signals from other cells and Z Un is the additive 
noise or extra interference. Z Un is modeled for convenience as white Gaussian with power 
E|Z U J 2 = a 2 . Note that a single power level is applied at each AP in this notation. This will 
allow us to ease the exposition of our analysis. In the OFDMA case however, a possibly unequal 
power level may be applied on each subcarrier, leading to the optimization of a power vector, 
under sum power constraint, rather than a scalar power level at each AP. The analysis in that 
case however leads to similar conclusions on the capacity scaling and is skipped in this paper. 



As stated above, intra-cell multiple access is orthogonal, while intercell multiple access is 
simply superposed, due to full reuse of spectrum. The resource allocation problem considered 
here consists in power allocation and user scheduling subproblems. Importantly we focus on 
rate maximizing resource allocation policies, rather than fairness-oriented ones. As is the case 
with known single cell protocols, multicell scheduling protocols can be enhanced to offer some 
desired performance-fairness trade-off, however this is outside the focus of this paper. Fairness 
issues are touched upon in [11]. In our setting the optimization of resource in the various resource 
slots decouples and we can consider the power allocation and user scheduling maximizing the 
capacity in any one slot, independently of other slots. A few useful definitions follow. 



N 




III. The multicell resource allocation problem 
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Definition 1: A scheduling vector U contains the set of users simultaneously scheduled across 
all N cells in the same slot: 

U = [u l U 2 ■ ■ ■ U n ■ ■ ■ U N ] 

where [U] n = u n . Noting that 1 < u n < U n , the constraint set of scheduling vectors is given by 
T = {U | 1 < u n < U n V n = 1, . . . , TV}. 

Definition 2: A transmit power vector P contains the transmit power values used by each 
AP to communicate with its respective user: 

P — [Pui P-U2 ' ' ' Pu n ' ' ' Pu N ] 

where [P] n = P Un = E|X M J 2 . Due to the peak power constraint < P Un < P max5 the constraint 
set of transmit power vectors is given by Vl = {P | < P Un < P max V n = 1, . . . , iV}. 

A. Capacity optimal resource allocation 

The merit (or utility) associated with a particular choice of a scheduling vector and power 
allocation vector is measured via the set of SINRs observed by all scheduled users simultaneously. 
T([U] n ,P) refers to the SINR experienced by the receiver u n in cell n as a result of power 
allocation in all cells, and is given by: 

T([U] n ,P) = , (1) 

°" 2 + ^2 Gu„,iPui 

where G Unti = \ a Unti \ 2 is the channel power gain from cell % to receiver u n . 

In data-centric applications, a reasonable choice of utility is a monotonically piece-wise 
increasing function of the SINR, reflecting the various coding rates implemented in the system. 
With an idealized link adaptation protocol, the utility eventually converges to a smooth function 
reflecting the user's instantaneous rate in Bits/Sec/Hz. For the overall network utility, and under 
the assumption that the various cell transmitters cannot afford to perform cooperative coding, we 
consider the average of rates achieved over all cells under single-user decoding, thus absorbing 
all sources of interference under a global Gaussian interference process. We may write the utility 
as [14]: 

1 N 

c(u,p) ^ - J2^g(i + r([u]n,P)). (2) 

71=1 
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The capacity optimal resource allocation problem can now be formalized simply as: 

(C/*,P*) = a rgmaxC(C7,P), (3) 
Pen 

The optimization above can be seen as generalizing known approaches in two ways. First the 
capacity maximizing scheduling problem has been considered (e.g. [15]), but in general not 
jointly over multiple cells. Second, the problem above extends the classical multicell power 
control problem (which usually rather aims at achieving SINR balancing) to include joint 
optimization with the scheduler. 

The problem in © presents the system designer with many degrees of freedom to boost 
system capacity but also with several serious challenges. First the problem above is non convex 
and standard optimization techniques do not apply directly. On the other hand an exhaustive 
search of the (U, P) pairs over the constraint set is prohibitive. Finally, even if computational 
issues were to be resolved, the optimal solution still requires a central controller updated with 
instantaneous inter-cell channel gains which would create acute signaling overhead issues in 
practice. The central question addressed by this paper can be formulated as follows: Can we 
extract all/some of the gain related to multicell resource allocation (compared with single cell 
treatment) within reasonable complexity and signaling constraints? Inspection of the recent 
literature reveals that this is a hot research issue with many possible tracks of investigation 
including use of modified capacity metrics, game theoretic approaches, reuse partitioning, power 
shaping and power quantizing (see e.g. [11] and references therein). Below, we investigate 
theoretical answers on this question by means of so-called scaling laws of the capacity, obtained 
via extreme value theory. This study reveals both surprising and promising answers. 

Interestingly, other work exists on analyzing the scaling law of capacity in interference-limited 
networks, including recently submitted [16]. In such work, a similar metric is used related to 
the sum capacity of simultaneously active links. Importantly though, in their work the scaling 
is in terms of growing number of links (or cells for a cellular network), rather than growing 
number of users per cell in a multiple access scheme as is studied here. Thus multiuser diversity 
scheduling is not exploited and fundamentally different scaling laws are obtained in the two 
cases. 
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IV. Network capacity: Models and bounds 

Let us consider a system with a large number of users in each cell. For the sake of exposition 
we shall assume U n = U for all n, where U is asymptotically large, while N remains fixed. 
We expect a growth of the sum capacity C(U*,P*) with U thanks to the multicell multiuser 
diversity gairo Thus we are interested in how the expected sum capacity scales with U. To 
this end we shall use several interpretable bounding arguments. We consider two channel gain 
models. The first considers a symmetric distribution of gains to all users from their serving AP. 
In the other one, an additional random distance-dependent path loss is accounted for. 

A. Bounds on network capacity 

The simple bounds below hold in the asymptotic and non asymptotic regimes as well. 
Upper bound: An upper bound (ub) on the capacity for a given resource allocation vector 
(not necessarily an optimal one) is obtained by simply ignoring intercell interference effects: 



i * , n p x 

c(tf,p)<^£iog(i + ^^). 



N 

i (,.. ./-'., \ 

(4) 



N " v a 

n=l 

In the absence of interference, the optimal capacity is clearly reached by transmitting at a level 
equal to the power constraint, i.e. ~P max = [P ma x, ■ ■ ■ , Pmax) and selecting the user with the 
largest channel gain in each cell (maximum rate scheduler), thus giving the following upper 
bound on capacity: 



C(U*,P*) < C ub (5) 



where 



1 N 

r " = ^E lo g( 1+r « 6 )- (6) 

71=1 

and where the upper bound on SINR Tjf is given by the maximum rate scheduler: 

Tl b = max {G Unin }P max /a 2 (7) 

u n =l...U 

'The multicell multiuser diversity gain can be seen as a generalization of the conventional multiuser diversity [15] to multicell 
scenarios with joint scheduling 



DRAFT 



9 



Lower bound: A lower bound (lb) on the optimal capacity (in the presence of interference) 
C(U* , P*) can be derived by restricting the domain of optimization. Namely, by restricting the 
power allocation vector to full power P ma x in all transmitters, we have 

C(U*,P*) > C lb (8) 

where 



C — C(U ppi Pmax) (9) 

and where U* FP denotes the maximum rate scheduling vector when assuming full power 
everywhere. This scheduling vector is defined by 



U* FP = argmaxC(£7, P max ), (10) 
Note that the user selected in the n-th cell, designated by [U* FP ] n , is found via: 

\TT* 1 {G Un ,n}Pmax f-\-\\ 

\U pp n = argmax >? (11) 

L FP\n 6 N 

u ~ / j\-f-r\ max 

The SINR corresponding to the selected user, denoted by T l £, is therefore given by: 



r'°= max L Z m " x (12) 

u n =l... U(T 2 + J2^ n G Un ,P max 

Finally the lower bound on capacity C lb may be rewritten as: 



N 

Mb 

L N 

n=l 



1 N 

i^iog(i + r'„'). (13) 



B. Distributed vs. centralized scheduling 

For large networks, it is important that scheduling algorithms can operate on a distributed 
mode, that is, the choice of the optimal user set should be done by each cell on the basis of 
locally available information only. This is in principle difficult task because the achievable rates 
observed in different cells are coupled together through the interference terms. However, we note 
that if the scheduler is based on maximizing the upper bound of network capacity given by ©, 
then each cell only needs to know the realization of the direct gain G Un>n , and the scheduler is 
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trivially distributed. Alternatively, to obtain a scheduler maximizing the lower bound of capacity 
given by (fT3l) . each cell must collect the worst case SINR for each of its users. The worst case 
SINRs are computed during e.g. a common preamble phase where all APs are asked to transmit 
pilot or data symbols at full power. This makes the scheduler of (fTTI) also distributed. Note that 
"worst case" is here understood in terms relative to the power control policy, not the scheduler. 

C. Channel models 

We now detail our assumptions regarding the fading and path loss models. Some of these 
assumptions are mainly technical, serving to simplify the analysis but could be relaxed without 
altering the fundamental results, as discussed later. As mentioned above we assume a cellular 
network where APs are regularly located with cell radius R. In this sense, the cells are assumed 
to be circular with each base being at the center of it, although this assumption is not critical to 
this study (i.e. similar conclusions can be obtained for hexagonal cell etc.) as explained below. 

The basic channel model consists in the product between a variable representing the path loss 
and a variable representing the fast fading coefficient: Let G Unji = r y Un ,i\h Un ,i\ 2 , u n = 1 . . .U,i = 
1 ... AT be the set of power gains where j Unji is the path loss between user u n (selected in cell n) 
and the access point in cell i. h Unji is the corresponding normalized complex fading coefficient. 
A generic path loss model is given by [17]: 

lu n ,i = PdZ,i (14) 

where (3 is scaling factor, e is the path loss exponent (usually with e > 2), and d Unji is the 
distance between user u n and AP i. 

Note that we assume as preamble a user-to-AP assignment strategy resulting in all users 
being served by the AP with the smallest path loss. This means, as is usually the case in current 
network design, that the AP assignment operates on a time scale which is not fast enough to 
provide diversity against fast fading. 

We consider in turn two basic user location scenarios, and a hybrid one. As it will be made 
clear later, the user location scenario has significant impact on the analysis of the network 
capacity. In the first scenario, denoted as symmetric network, all users served by a given AP 
are assumed to be located at the same distance from that AP. This idealized situation results in 
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all users experiencing the same average SNR, an assumption often made by previous authors in 
this area, and for which several interesting results of the existing literature can be reused. This 
scenario is illustrated in Figj2l 

In the second, more realistic, scenario, denoted simply as non-symmetric network, the users 
are located randomly over a cell given by a disk of radius R around each of the serving APs. 
Finally, a hybrid scenario mixing the two scenarios above is discussed later in the paper. 

Note that the actual cell shape will not be a disk in reality. However we argue that, when it 
comes to studying the scaling laws of network capacity with maximum-capacity user scheduling, 
the actual shape taken by the cell borders has in fact little impact on the result. The main reason 
is that since the user's direct links is subject to a location dependent path loss, the distance to 
the serving AP will affect its chances of being selected by the scheduler. As a consequence the 
users located in the inner region of the cell (i.e. close to the access point) bear the vast majority 
of the traffic and are the drivers for the capacity scaling laws. Therefore an accurate modeling 
for the location of cell-edge users is unimportant here. 

V. Network capacity: Scaling laws 
A. Capacity scaling with large U in symmetric network 

We analyze the scaling of capacity C(U*, P*) via the scaling of the bounds C lb and C ub , with 
increasing U. We just focus on the performance in cell n, as other cells are expected to behave 
similarly under equal number of users U and isotropic conditions throughout the network. For 
the symmetric network, users experience an equal average SNR, thus j Un , n = In is a constant 
independent of the user index. 

Interestingly, for this particular case, we show we can reuse extreme value theory results [18] 
developed specifically in the context of single cell opportunistic beamforming [12], [13] and 
transposed here to the case of networks with multicell interference. For the case asymmetric 
network, specific results are developed. 

First, the following results provide insight into the interference-free scaling of SINR and 
capacity respectively. 

1 ) Scaling laws for interference-free case: The interference-free multicell capacity scaling 
boils down to studying the scaling in each cell independently. Further assuming an isotropic 
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network (i.e. all cells experience the same channel statistics) we can simplify the analysis by 
exploiting known results on single-cell capacity scaling, as done below. 

Lemma 1: Let G UniJl = 7«n,nl^w| 2 > u n = 1 . . . U, n = 1 . . . N, where j Unjn = j n . Assume 
\h Un ,n\ 2 is Chi-square distributed with 2 degrees of freedom (% 2 (2)) (i.e. h Un>n is a unit-variance 
complex normal random variable). Assume the \h UnyTl \ 2 are independent and identically distributed 
(i.i.d.) across users. Then for fixed N and U asymptotically large, the upper bound on the SINR 
in cell n scales like 

Vf « *™LhL]ogU (15) 

where the symbol ps means that the ratio of the left hand side and right hand side terms converges 
to one, as U goes to infinity. 

Proof: This result is a reuse of a now well known result for single cell opportunistic scheduling. 
This states that the maximum of U i.i.d. x 2 (2) random variables behaves like log £7 for large U. 
See for instance [12], itself building on classical extreme value theory results [18]. We omit the 
proof here and refer the readers to these references. 

From the scaling the SNR, we obtain the scaling of the interference-free capacity shown in 
|6] This is stated in the following theorem, again building on known single cell results but stated 
here for convenience, with our own notations: 

Theorem 1: Let G Un>n = 7wl^wl 2 > u n = l...U,n = 1...N, where 7 W = 7. This 
means that all cells are assumed to enjoy an identical link budget. Assume \h Untn \ 2 is Chi-square 
distributed with 2 degrees of freedom (x 2 (2)). Assume the \h Un ^ n \ 2 are i.i.d. across users. Then 
for fixed N and U asymptotically large, the average of the upper bound on the network capacity 
scales like 

E(C ub ) w log log [/ (16) 

where the expectation is taken over the complex fading gains. 
Proof: Under isotropic network conditions, we have from ©: 

E(C ub ) =E(log(l + r^)) (17) 

Once the scaling of Tjf is obtained, the scaling of the expected value of log(l + ) is readily 
obtained from published results in the context of single-cell maximum rate user scheduling, found 
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in [12], [13] among others. For a detailed proof, see e.g. [13], Theorem 1. 

2) Scaling laws for full-powered interference case: We now turn to the behavior of interfer- 
ence limited networks by exploring the lower bounds given for SINR and capacity. The initial 
intuition would be that the analysis of the lower bound given in (PT21) provides us with a very 
pessimistic view of the network performance as it assumes interference coming at full power 
from every AP in the network. The interesting aspect behind our findings below is that it is 
not. In fact the negative impact of interference at the user on network capacity can be made 
arbitrarily small while not sacrificing transmission rates to the assigned APs, as shown per the 
following theorems. In the results below, remember we assume each user is assigned to a serving 
AP which is the one with minimum path loss. As a consequence, since the region of coverage 
under study is limited to a disk of radius R around the serving AP, the distance between a user 
and any interfering AP is greater than R. As a result we have from (fl4l) : 

G Unti < (3R~ e \K r J 2 foranyi^n (18) 
The lemma below gives the scaling law for the worst case SINR [T2l 

Lemma 2: Let G Un)i = 'j Un ,i\h Unji \ 2 , u n = l...U,n = 1...N, where j Un;U = 7„, j Un;i = 
(3d~ e { for i ^ n. Assume \h Unji \ 2 is Chi-square distributed with 2 degrees of freedom (% 2 (2)). 
Assume the \h Unyi \ 2 are i.i.d. across users, cells. Then for fixed N and U asymptotically large, 
the lower bound on the SINR in cell n scales like 

-nib Pmax^n 



a 2 



\ogU (19) 



Proof: To obtain this result, one uses the fact that users in cell n are served by their closest AP. 
Following (TT8l . an upper bound on the interference power then given by YuiLn l3R~ e \h Un ^\ 2 P max . 
This gives a further lower bound on Y 1 ^ given by 

rt > rf (20) 

where T l £ 2 is corresponds to the SINR assuming pessimistically that all sources of interferences 
are located on the edge of the cell of interest, calculated by: 



^62 



"fnPmax ™aX U Un (21) 
u„=l...U 
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where u Un denotes the normalized SINR at user u 

\h 



uj u = r "" n ' n| „ (22) 

The scaling law of T^ 2 is also that of tu Un , which is the ratio of a Chi-square (2 degrees of 
freedom) distributed variable and the sum of a fixed noise term and a Chi-square (2N-2 degrees 
of freedom) variable. Thus the scaling of uj Un is similar to the scaling of the SINR in the single 
cell opportunistic beamforming problem with N antennas at the transmitter, studied in [13]. In 
there, the SINR is the ratio of a direct beam power and a noise plus N — 1 interfering beam 
power term. In particular we can find its distribution as: 



Fw{u}) ~ 1 " (lTW^Fp (23) 

([13], Lemma 4) shows that the SINR then scales like \ogU . This gives in our context: 

rf ^P maxln \ogU/a 2 (24) 

Note that the scaling above is identical to the one reported for the interference-free case (fl"5l) . 

Thus, is bounded above and below by two expressions (respectively the interference-free 
T^f and T^ 2 ) which exhibit the same scaling law. Therefore must satisfy itself the same 
scaling law. 

The following theorem gives the scaling law for the lower bound on capacity for an isotropic 
network. 

Theorem 2: Let G Un>i = j Unji \ h Unti | 2 , u n = l...U,n = 1...N, where 7 W = 7, j Unti = 
/3d~^ { for i 7^ n. Assume \h Un:i \ 2 is Chi-square distributed with 2 degrees of freedom (% 2 (2)). 
Assume the \h Un! i\ 2 are i.i.d. across users, cells. Then for fixed and U asymptotically large, 
the average of the lower bound on the network capacity scales like 

E(C lb ) w log log [/ (25) 
Proof: From the result in Lemma [2l this result is proved in a way identical with that of ([13], 
Theorem 1). Therefore the proof is omitted here for space considerations. 

From bounding arguments and from theorems [I] and [2] above, the following conclusion is now 
obtained: 
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Theorem 3: Let G UnA = j Un>i \ h UnA \ 2 , u n = 1 . . . U, n = 1 . . . N, where 7 Un , n = j n , 7 Unji = 
/3rf^i for i ^ n. Assume \h Unji \ 2 is Chi-square distributed with 2 degrees of freedom (% 2 (2)). 
Assume the |/i,i n) j| 2 are i.i.d. across users, cells. Then for fixed N and U asymptotically large, 
the average of the network capacity with optimum power control and scheduling scales like 

E(C(U*,P*)) « log log U (26) 
Proof: The result is readily obtained from writing: 

E(C lb ) < E(C(U*, P*)) < E{C ub ) (27) 

Then, invoking (T2~5l) and (fl6l) exhibiting the same scaling law, we obtain a similar law in (|26|) . 

Theorems \T\ and [2] suggest that, in a multicell network with symmetric users, the capacity 
obtained with optimal multicell scheduling in both an interference-free environment and an 
environment with full interference power have identical scaling laws in log log U. This result 
bears analogy to the results by [13] which indicate that in a single cell broadcast channel 
with random beamforming and opportunistic scheduling, the degradation caused by inter-beam 
interference tends becomes negligible when the number of users to choose from becomes large. 
Here the multicell interference becomes negligible because the optimum scheduler tends to 
select users on an instantaneous basis who have both a strong direct link to their serving AP and 
small interfering links from surrounding APs. Interestingly, the minimization of the multicell 
interference term should take away some degrees of freedom in choosing the users with best 
direct links, however not sufficiently so to affect the overall capacity scaling. 

Another interpretation of this result is in terms of our ability to find distributed scheduling 
schemes for maximizing the network capacity. The optimal multicell scheduler and power control 
solution would be hard to implement in practice. However from the observations above, a simple 
scheme based on each cell measuring the worst case SINR of each of its users (during e.g. a 
preamble) and selecting the users with the best worst case SINR as per (fT2~l) . will result in an 
quasi optimal behavior asymptotically (again, from a scaling perspective). Such a scheme does 
not require any exchange of information between the cells and the worst case SINR can be 
measured in one shot by each user and fedback to its serving AP. 

These results come as a complement to previously reported findings [19], [16] which propose 
a near optimal power allocation scheme, for fixed number of users, where a fraction of the 
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transmitters are selected to be turned off while the rest operate at full power. It was observed 
experimentally [19] there that the fraction of off cells would go to zero when the number of 
users grows large. Thus in a network with full reuse and greedy user scheduling, the optimal 
power control policy should be for all cells to operate at the power constraint. The analysis of 
scaling of capacity provides a theoretical justification to this intuitive result. 

We now turn to a non symmetric network where users can experience different average SNR 
values depending on their position and conduct a similar analysis. However we will see that 
different capacity scaling rates are obtained compared with the symmetric network case. 

B. Capacity scaling with large U in non-symmetric network 

We assume the path loss is determined by the user's distance to the emitting AP, both serving 
and interfering. We consider a uniform distribution of the population in each cell. Thus d Un>n 
(distance between user u n and its serving AP) is a random variable with non-uniform distribution 
For a cell radius R, we find easily: 

f D (d) = 2d/R 2 , de[0,R] (28) 

Further, the random process d Unyn can be considered i.i.d. across users and cells, if users in each 
cell are dropped randomly in each dislo Assuming R = 1 for normalization, the distribution of 
7«„,n = Pd~^i is given by (details omitted here): 

fM = { eP 9 (29) 
I with g £ [(3, oo) 

In order to get upper and lower bounds on performance, we are interested in the behavior of 
the following extreme values of product of independent random variables: 

max 1u n ,n\h Un , n \ 2 for the interference-free case and 

u n =l...U 

max ^u n ,n^u n for the full-powered interference case 

«n=l— U 



where u Un is again defined as per (|22|) . 

2 The considered coverage region can be assimilated with the inside area of each disk, in a disk-packing region of the 2D 
plane. Users dropped outside the disks can dropped from the analysis, as these will not affect the scaling law. 
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1 ) Extreme values of heavy -tail random variables: The distribution of 7 Mnj „ shown in (|29l ) 
is remarkable in that it differs strongly from fast fading distributions, due to its heavy tail 
behavior. Tail behavior clearly plays a fundamental role in shaping the limiting distribution of 
the maximum value, hence also the scaling of capacity. Note that heavy tail is also observed 
in large scale fading models such as log normal shadowing for instance. In order to study the 
extreme value of a product of random variables involving one heavy tailed variable, we need 
first to review the properties of so-called regularly varying random variables. See e.g. [18] for 
a definition of such variables, restated below: 

Definition 3: A random variable X, with distribution (cdf) given by F x (x), is said to be 
regularly varying (at oo) with exponent —a if and only if: 

1 - F x (x) ^ when _^ 
l-F x (tx) 

The lemma below shows how the definition above applies to our situation: 

Lemma 3: Let X = 7« njn . X is regularly varying with exponent — -. 

Proof: A direct application of the definition above, with a distribution obtained from (T29l) : 

F x (x) = 1 - Q' 2 - x>p. (31) 

An interesting aspect of regularly varying distributed random variable (R.V.) is that they are 
stable with respect to multiplication with other independent R.V. with finite moments as pointed 
out by the following theorem shown by Breiman [20]: 

Theorem 4: Let X and Y be two independent R.V. such that X is regularly varying with 
exponent —a. Assuming Y has finite moment E(F a ), then the tail behavior of the product 
Z = XY is governed by: 

l-Fz(z) ^E(Y a )(l- F x (z)) whenz^oo (32) 
The idea behind this theorem is that when multiplying a regularly varying R.V. with another 
one with finite moment, one obtains a heavy tailed R.V. whose tail behavior is similar to the first 
one, up to a scaling. In other words, heavy tail behavior tends to dominate over other distribution. 

We now apply this result to X = 7« n , n and Y given by Y = \h Unin \ 2 for the interference free 
case and Y = cu Un for the full-powered interference case, respectively. Note that in both cases, 
Y has finite moments. The tail behavior of Z — XY can then be characterized by the following 
lemma: 
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Lemma 4: Let X = j Unjn be a R.V. with distribution given by (|29l ). Let Y be an independent 
R.V. such that E(Y?) < oo. Then the tail of Z = XY is governed by: 

1 - F z (z) -> E(y «) whenz^oo (33) 

Proof: A direct application of Theorem |4] using the distribution of X shown in |3~T1 

The lemma above indicates that the tail behavior of the distribution of X = 7 11n , n , characterized 
by Lemma [3] carries over to that of the product Z = XY. As a consequence, Z is also regularly 
varying with the same exponent — |. 

We now complete our study by reviewing existing results on the extreme value of regularly 
varying R.V. Following [18], a regularly varying variable can be classified to be of Frechet 
type. Extreme values of Frechet (or regularly varying) variables are characterized by use of 
the Gnedenko theorem, given in appendix I. For comparison, note that the random variables 
involved in the analysis of previous sections (Sec lV-Al and therein), belong to the so-called 
Gumbel category. In our context, we have the following result: 

Lemma 5: Let Z Un = 7 nni „F where Y is a R.V. with finite moments, independent of 7 Unjn . 
Then we have: 

limPr{ max Z Un < /M(Y*)%Uh} = e _i " ? Vt > 0, when U -> oo (34) 

u n =l...U 

Proof: We invoke Gnedenko's theorem [21] given in appendix I. It is easy to find that au = 

/?E(Kv)2[/2 where ajj is defined in the appendix. 

2) Scaling law for interference-free case: The inequality in © allows us to characterize the 
scaling law of capacity. Although a characterization in terms similar to those of previous section 
(i.e. finding a scaling law l(U) for the SINR, such that the ratio of the SINR and l(U) converges 
towards 1 when U — > oo) is possible when analyzing the capacity, such a task is very tedious and 
mathematically involved. For the sake of easing exposition, we somewhat loosen our definition of 
scaling of SINR below in a way that allows us to directly exploit Lemma \5\ while preserving the 
key interpretations. Briefly speaking, we characterize the scaling of SINR up to a multiplicative 
constant. However it's important to note that the scaling of capacity itself, within the standard 
definition of scaling, will be determined precisely (thanks to the log function behavior). 

The theorem below gives the scaling law of SINR for the interference-free case in a non- 
symmetric network. First we give the following wider- sense definition of scaling: 
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Definition 4: Let g{U) and l(U) two functions, defined over U G [C/ ,+oo], where U is an 
arbitrary real. g(U) is said to scale as l(U), i.e. g(U) ~ 1{U), U — > oo when 

Pr(g(U) > v(U)) -> 0, when [/ -»• oo 

Pr(#(£f) < tu(Ef)) -> 0, when ?7 -»• oo 

(35) 

for any two functions v(U) and w(U) such that — > and ^© — > 0, respectively. 

Note that this notion of scaling can be interpreted as g(U) grows neither significantly faster 
than l(U), not does it grow significantly slower than l(U). 

Theorem 5: Let h UnjU , u n = 1 . . . U be i.i.d. Gaussian distributed unit-variance random vari- 
ables. Assuming that 7« njn is i.i.d., distributed as per (|29l) , for n = 1 . . . N. Then for fixed iV 
and U asymptotically large, the interference-free SNR scales like: 

K b ~ C/f (36) 
Proof: Let v(U) be any function growing faster than U^, i.e. such that lim^oo U%/v(U) = 0. 
Then let t = v(U)/(f3E(Y-)^U% ). From Lemma [5] we have that 

Pr{ max Z Un < v(U)} -> lim e"^ 1 = 1 (37) 

Equivalently, we have that Pr(max Uri= i...[/ Z Un > v(U)) — > 0. Similarly, we can prove that any 
function u>([7) growing slower than U% will be such that Pr{max Un=1 ><<t / Z Un < w (£/")} — > 0. 
Thus max Un=1 j/ Z Un scales as f/i in the sense of definition 4. 

From the scaling of SNR above, we can infer the exact scaling on the upper bound on capacity 
E(C ub ), as shown per the theorem below: 

Theorem 6: Let h UnjU , u n — 1 . . . U be i.i.d. Gaussian distributed unit- variance random vari- 
ables. Assuming that j Un ,n is i.i.d., distributed as per (|29l) , for n = 1 . . . N. Then for fixed iV 
and U asymptotically large, the interference-free capacity scales like (i.e. the ratio of the two 
quantities converges to 1 almost surely): 

E(C ub ) w -logt/ for large U (38) 

Proof: See appendix II. 

We now proceed to determine the scaling laws in the case of full-powered interference. 
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3) Scaling law for full-powered interference case: We can derive the scaling laws for the 
lower bound of SINR and capacity by following a strategy similar to Sec lV-B.2l simply by 
replacing the R.V. \h Un ^ n \ 2 by the R.V. uj Un which also has bounded moments. We obtain the 
following result: 

Theorem 7: Let h Unti , u n = 1 . . . U, i = 1 ... N be i.i.d. Gaussian distributed unit-variance 
random variables. Assuming that 7 Uni „ is i.i.d., distributed as per (l29l) . for n — 1 . . . N. Then for 
fixed iV and U asymptotically large, the lower bound on SINR scales like: 

It ~ Ui (39) 
Proof: We use the same proof as for Theorem [51 with X = 7 Uni „ but this time Y = cu Un . 

Finally, from Theorem [71 we infer that the upper bound on capacity for a non-symmetric 
network exhibits an exact scaling defined as: 

E(C lb )^UogU (40) 

The proof for ([401) is identical to that of Theorem [6] in Appendix II, but simply replacing Y 
with cu Un , which clearly does not change the scaling. 

Remarkably, as in the case of the symmetric network, the results above (l38l) and (|40"1) suggest 
that multicell interference, no matter how strong, does not affect the scaling of the network 
capacity, if enough users exist and rate-optimal scheduling is applied. Furthermore, by virtue of 
the upper bound and lower bound exhibiting the same scaling law in (l38l) and (l40l) respectively, 
the capacity under optimal scheduling and power allocation must behave like 

c(tr,p*) «| log c/ (4i) 

Two remarks are in order. First, in the symmetric network case, a suboptimal but fully 
distributed resource allocation based on constant (full) power transmission at all transmitters 
and scheduling policy based on (TTTI) will actually result in the best possible scaling law of 
capacity for the network. Second, we observe that we obtain a much greater rate growth than 
in the case of the symmetric network. This is due to the amplified multiuser diversity gain due 
to the presence of unequal path loss across the user locations in the cell. This results from a 

DRAFT 



21 



scheduler which, in a quite unfair fashion admittedly, tends to select users closer to the access 
point as more users are added to the network. 

C. Discussion on channel models and exclusion area around the AP 

Interestingly, the theory on regularly varying variables stipulates that multiplication of the path 
loss variables by any small scale fading variable with finite moments will preserve its heavy tail 
behavior. This means that our result shown in (|4Tj) is in fact valid for a wider class of fading 
channel models, such as Nakagami, Rice, etc. On a different note, one may wonder how close 
users can be assumed to get to the access point in practice. Let us imagine that a small disk of 
exclusion, with the AP at its center, prevents users to getting to close to the AP. As a by product, 
the disk also serves the purpose of maintaining the validity of the path loss model, which may 
not be reasonable in the close vicinity of the AP. In this case, one may expect two successive 
regimes for the capacity scaling as U grows. In the first regime, when the number of users is 
still moderate, the scaling is dominated by the path loss effect, with a behavior such as shown 
in (|4TT) . In the second regime, when enough users are already accumulated near the exclusion 
circle, it is the turn of the tail behavior of small scale fading to dominate and the scaling will 
be characterized by (|26l) . This situation is investigated in the next section. 

As the growth would be ultimately limited by that the tail of the small-scale random fading in 
practical situations, one may also wonder how accurately Chi- square distributions model reality 
in real-world wireless channels. Clearly, this discussion is inherent in all previous studies dealing 
with scaling laws and asymptotic performance analysis. Nevertheless it is important to keep in 
mind the basic law of power preservation which indicates that no matter how many users are 
considered, the most favorable users can't receive more power than what was actually transmitted. 
This simple fact will impose a hard limit on the SNR which in turn limits the domain of validity 
of our scaling in terms of the number of users U. Although we believe a specific analysis of the 
validity domain will rely on yet unexpored channel model properties (tail properties of the pdf 
are less explored than the behavior near zero which characterize outage) and is outside the scope 
of this paper, it remains clear that this domain is wide enough for the analysis to be meaningful 
since the power preservation limit is reached only when the small scale fading is in the order of 
the inverse of path loss, which would require very large fading coefficients in practice (several 
tens of dB). 
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VI. Numerical evaluation 

We validate the asymptotic behavior of the multicell sum rate when U grows large with 
Monte Carlo simulations. We use a network with N = 4 cells, unit cell radius and the following 
parameters @ — 1/16, e = 4, P max = 1, a 2 = 0.02. I.i.d. flat Rayleigh fading is considered 
in addition to the path loss based power decay. We consider three scenarios for user location, 
as mentioned previously in this paper. First, we consider cells with users located on a circle 
with distance 0.5 away from the AP (symmetric network). Then we consider a non-symmetric 
distribution of average SNR by drawing users randomly in the cell. Finally we consider an hybrid 
scenario where users are drawn unform randomly over the cell but kept outside an exclusion 
disk of radius 0.1 around the AP. In all cases, we evaluate the upper and lower bound on per-cell 
data rates (see FigJH FigJH Figj5] and observe the identical rate growth of the lower and upper 
rate bounds. This also shows that the capacity obtained with exhaustive user and power level 
selection also has the same growth rate. The observed rate growth in log log U for the symmetric 
network and in log U for the non symmetric network confirms our earlier theoretical claims. In 
FigjH we observe a scaling behavior with two distinct regimes with a \ogU in the moderate 
number of users U and log log U for high number of users, thus confirming our intuition for 
what could happen in a realistic network. 

VII. Conclusions 

We present an extreme value theoretic analysis of network capacity for maximum sum rate 
multicell power allocation and user scheduling. We derive scaling laws of capacity when the 
number of users per cell grows large, both in cases where the users have same average SNR 
and path loss dependent SNR. We show that in both cases, 1-the effect of intercell interference 
on rate scaling tends to be negligible asymptotically, and 2- should intercell interference be 
considered, an asymptotically optimal allocation procedure is given based on full power allocation 
at all transmitters, which is furthermore completely distributed. We show that the growth rate 
of capacity is exponentially faster in the case of a system with unequal distance-based average 
SNR. 
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Appendix I 

The following theorem is due to Gnedenko [21] (1943) and states the following property for 
regularly varying distributions: 

Theorem 8: Let Zi an i.i.d. random process. Then Z^ has a regularly varying distribution with 
exponent a if and only if 

limPr{ max Z< < a v t} = e~ l ~ a Vt > when U -> oo (42) 

i=l. ..[/ 

where a;y is a sequence such that 1 — F z {a u ) = jj. 

Appendix II 

From Lemma [5] we have that 

limPr{r^ < (3E(Y~*)%Uh} = e"*" 1 Vt > 0, when U -> oo (43) 

Since the SNR rjf is growing large in each cell by virtue of Theorem \5\ the capacity can be 
approximated by: 

1 - 

n=l 

when [/ grows large. From (|43T) . we write 

limPr{logrf < bg(/3E(yt)f ) + log t + | logC/} = e"*" 1 Vt > 0, when [/ -»■ oo (45) 
Now, taking £ = log [/ we infer that 

logl^ 6 < log(/3E(y«)5) + loglog£/ + ^logf/ almost surely when C/ -> oo (46) 
On the other hand, taking t = 1/ log [/, we obtain that 

log rjf > ^ log U — log log [/} almost surely when U — > oo (47) 

From (1461) and (1471) . we conclude that 

log T ub 

> 1 almost surely when U — > oo (48) 

I log 17 

From the isotropy of the network, this shows that C ub (and a fortiori its average) scales as 
flogt/. 
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Fig. 1. A two-cell network diagram example. Direct and interfering links toward the scheduled user (black) are indicated in 
solid and dashed arrows respectively. Users are located randomly over a cell of radius R around their access point. 




Fig. 2. A two-cell idealized symmetric network diagram example. Direct and interfering links toward the scheduled user (black) 
are indicated in solid and dashed arrows respectively. In this idealized case, users a located over a circle, a fixed distance away 
from their access point. 
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3. Scaling of upper and lower bounds of capacity versus U for a symmetric network (N = 4). The observed scaling for 
curves is in log log U. 
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Fig. 4. Scaling of upper and lower bounds of capacity versus U for a non-symmetric network (unequal average SNR) (TV = 4). 
The observed scaling for both curves is in log U. 
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Fig. 5. Scaling of upper and lower bounds of capacity versus U for a non-symmetric network with a an exclusion disk around 
each AP of radius 0.1. (N — 4).The observed scaling in the transitory regime is in log U, then in log log U in the asymptotic 
regime. 
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